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Abstract 

We present results of the clustering analysis between active galactic nuclei (AGNs) and galax¬ 
ies at redshift 0.1-1.0, which was performed for investigating properties of galaxies associated 
with the AGNs and revealing the nature of fueling mechanism of supermassive black holes 
(SMBHs). We used 8059 SDSS AGNs/OSOs for which virial masses of individual SMBHs were 
measured, and divided them into four mass groups. Cross-correlation analysis was performed 
to reconfirm our previous result that cross-correlation length increases with SMBH mass MbhI 
we obtained consistent results. A linear bias of AGN for each mass group was measured as 
1.47 for Mbh = M© and 3.08 for Mbh = IO^-IO^^Mq. The averaged color and 

luminosity distributions of galaxies around the AGNs/QSOs were also derived for each mass 
group. The galaxy color Hopt-m was estimated from an SED constructed with a merged SDSS 
and UKIDSS catalog. The distributions of color and luminosity were derived by a subtraction 
method, which does not require redshift information of galaxies. The main results of this work 
are: (1) a linear bias increases by a factor two from the lower mass group to the highest mass 
group; (2) the environment around AGNs with the most massive SMBH (Mbh > 10® Mq) is 
dominated by red sequence galaxies; (3) marginal indication of decline in luminosity function 
at dimmer side of Mir > -19.5 is found for galaxies around AGNs with Mbh = 10®-^ - 10® M© 
and nearest redshift group (z =0.1-0.3). These results indicate that AGNs with the most mas¬ 
sive SMBHs reside in haloes where large fraction of galaxies have been transited to the red 
sequence. The accretion of hot halo gas as well as recycled gas from evolving stars can be 
one of the plausible mechanisms to fuel the SMBHs above ^ 10® Mq. 

Key words: astronomical databases: miscellaneous — galaxies: active — large-scale structure of uni¬ 
verse — quasars: general — virtual observatory tools 


1 Introduction 

There are a lot of observational evidences that a Supermassive 
Black Hole (hereafter SMBH) is located in the center of all but 


the smallest galaxies (Richstone et al. 1998). Although the evo¬ 
lution mechanism of SMBH is still not well known, a recent 
growing evidence suggests that there is a strong link between 
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the growth of SMBH and the star formation in the host galaxy. 
One of the observational evidence is the similarity between the 
evolutions of black hole accretion rate and the star formation 
rate of galaxies (e.g., Madau et al. 1996; Boyle et al. 1998; Ueda 
et al. 2003; Zheng et al. 2009). Another important evidence is 
the correlation between mass of the SMBH (Mbh) and mass 
(Mh) / velocity dispersion (a) of the bulge component of its 
host galaxy (Magorrian et al. 1998; Ferrarese & Merritt 2000; 
Gebhardt et al. 2000; Ho 2007). These relations indicate that 
BHs and bulges coevolve regulating each others growth (see 
the recent review by Kormendy & Ho 2013, and references 
therein). 

According to the theoretical works (e.g., King 2014), mass 
of a SMBH can be regulated by its own outflow from an active 
galactic nucleus (AGN) in a way that the interstellar gas is ex¬ 
pelled far away from the host galaxy once Mbh reaches to the 
critical mass given by the Mbh-ct relation. As a result both the 
star formation in the bulge and mass accretion onto the SMBH 
are terminated. Therefore understanding the evolution mecha¬ 
nism of SMBHs and AGN phenomenon is crucial to shed light 
on the evolution of galaxies. 

The most important open questions are the nature of the fu¬ 
eling of SMBHs and triggering mechanisms of AGNs. The 
following three modes have been proposed for transferring gas 
onto the center of galaxy (e.g., Croton et al. 2006; Lagos et al. 
2008; Keres et al. 2009; Fanidakis et al. 2013): 

First one is a secular mode which arise through internal dy¬ 
namical processes in the disk such as a bar instability or external 
processes of galaxy interactions. Some authors claim that low 
luminosity AGNs hosted by late type disk galaxies are driven 
by this mode (see review by Kormendy & Kennicutt 2004, and 
references therein). 

Second one is a merger mode in which gravitational torques 
induced by galaxy-galaxy major/minor mergers drive inflows 
of cold gas toward the center of galaxies, triggering the central 
starbursts and also accretion on to the SMBH, i.e. AGN (e.g., 
Hopkins et al. 2008). There are observational evidence that 
some of the luminous AGN, i.e. QSOs, are likely triggered by 
a major merger (e.g., Sanders et al. 1988; Treister et al. 2012). 

There has been an observational evidence that AGN activity 
is enhanced for close companion galaxies which are assumed 
to be undergoing early stage of interaction (Silverman et al. 
2011). However, the fraction of such close pairs is only ~18%, 
and visual inspection of HST image revealed that ~ 85% of 
AGN host galaxies show no strong distortions on their mor¬ 
phologies (Cisternas et al. 2011). Thus galaxy interaction in¬ 
cluding major merger cannot be the only mechanism for fueling 
SMBHs, at least for AGNs with moderate luminosity sampled 
by Silverman et al. (2011) and Cisternas et al. (2011). Kaviraj 
(2014) analyzed SDSS data to probe the role of minor mergers 
in driving stellar mass and BH growth in galaxies, and suggest 


that around half of the star formation activity is triggered by the 
minor-merger process. 

Third one is a hot halo mode which is characterized by qui¬ 
escent accretion of gas from the hot halo (Keres et al. 2009; 
Fanidakis et al. 2013). In this mode the accretion rate is lower 
than the other two modes. However it becomes increasingly 
effective accretion mode at higher dark matter halo masses of 
Mh > 10^^ ® Mq, as the accretion by the other two modes be¬ 
comes inefficient due to AGN feedback. Thus larger fraction 
of SMBHs are expected to evolve through the hot halo mode as 
the host halo mass increases above Mh > 10^^'® Mq, and the 
corresponding mass of SMBH is > 10® Mq (Fanidakis et al. 
2013b). 

Due to the lower accretion rate, the expected AGN lu¬ 
minosity in this mode is typically lower than those of the 
other two modes in a dark matter halo with the same mass. 
Observationally such accretion is often associated with BHs that 
are characterized by very low radiative efficiency. 

Kauffmann et al. (2009) observationally found that there are 
two distinct modes in BH growth. One is associated with galaxy 
bulges that are undergoing significant star formation, and the 
other is associated with those with little or no on-going star for¬ 
mation. The former one could be related with the star formation 
induced by secular and/or merger mode, and the latter with the 
hot halo mode. 

In addition to the three modes, a model of accretion of recy¬ 
cled gas from evolving star has been studied (Ciotti & Ostriker 
2001, 2007). According to this model, the recycled gas accu¬ 
mulates within a galaxy under the influence of gas heating from 
Type la supernovae, and a fraction of the gas intermittently ac¬ 
cretes to the central SMBH under the effect of AGN feedback. 
With this scenario the central SMBH grows up to 10®'® Mq in 
giant ellipticals. 

The mass of dark matter halo in which SMBH resides can be 
estimated through the large-scale environments. As an estima¬ 
tor of large-scale environments, cross-correlation length and/or 
a bias parameter are usually used. The absolute bias of AGNs 
distribution relative to that of the dark matter can be derived 
from the observed two point correlation function, which is then 
used to estimate mass of the dark matter halo. 

Ross et al. (2009) performed correlation analysis for SDSS 
QSOs with redshift from 0.3 to 2.2, and obtained a result 
that the QSOs inhabit dark matter haloes of constant mass of 
~ 2 X 10^® Mq at any redshifts. Combining this observa¬ 
tional result with the theoretical work by Fanidakis et al. (2013), 
the optical AGNs/QSOs are expected to be mostly powered by 
either secular or merger mode. It is also inferred that they ap¬ 
pear only in the dark matter halo of which the mass does not 
exceed a critical mass that is determined by the AGN feedback. 

Hickox et al. (2009) analyzed the data of radio. X-ray, and 
IR detected AGN samples, and found that the correspond- 
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ing masses of dark matter haloes are ~ Mq, ~ 

Mq, and ~ Mq, respectively. Krumpe 

et al. (2012) also derived the halo mass for the X-ray selected 
ROSAT AGNs and the optically selected SDSS AGNs, and they 
obtained ~ Mq and ~ Mq, respectively. 

Shen et al. (2011) studied the dependence of quasar clus¬ 
tering on luminosity, mass of SMBH, quasar color, and radio 
loudness. They did not find significant dependence on either lu¬ 
minosity, mass, or color, while they marginally found that the 
most luminous and most massive QSOs are more strongly clus¬ 
tered than the remainder of the sample at ~2 sigma level, and 
radio-loud QSOs are more strongly clustered than radio-quiet 
QSOs. 

Many other studies on clustering and/or environment of 
AGNs/QSOs have been reported elsewhere (e.g.. Groom et al. 
2005; Coil et al. 2009; Silverman et al. 2009; Donoso et al. 
2010; Allevato et al. 2011; Bradshaw et al. 2011; Mountrichas 
et al. 2013; Zhang et al. 2013; Georgakakis et al. 2014). 

All these studies indicates that some types of AGNs such 
as those bright in the radio and/or X-ray band inhabits rela¬ 
tively larger dark matter haloes than optically bright AGNs, and 
at least part of them may be fueled by a different mechanism 
such as hot halo mode or accretion of recycled gas from evolved 
stars. 

To investigate the relation between the environment and in¬ 
trinsic properties of AGNs with better accuracy, we have devel¬ 
oped a clustering analysis method which does not require mea¬ 
surements of galaxy redshift and thus can utilize all the galaxies 
detected in the imaging data (Shirasaki et al. 2011; Komiya et 
al. 2013). Komiya et al. (2013) examined dependence of AGN- 
galaxy clustering on BH mass, and found an indication of an 
increasing trend of cross-correlation length above 10®’^ Mq. 

In this paper, following the result of Komiya et al. (2013), 
we study the properties of galaxies, such as color and lumi¬ 
nosity function, around AGNs and its dependence on BH mass 
by combining optical and near infrared dataset of SDSS and 
UKIDSS survey. We also derived AGN bias parameter to be 
compared with other studies. 

It is expected that, if the increase of clustering strength for 
AGNs with higher BH mass is due to the prominence of hot 
halo mode, most of the galaxies around the AGNs are in red 
sequence since they reside in haloes where gas cooling and star 
formation have been shut off by AGN feedback. Accretion of 
recycled gas from evolving stars may overcome in such a cir¬ 
cumstance. Which accretion occurs may depends on the avail¬ 
ability of the gas accumulated inside the galaxy and the strength 
of the AGN feedback. As another possibility, if the AGNs with 
the most massive SMBH are mostly fueled by mergers with gas- 
rich galaxies, the environment of the AGNs may be dominated 
with gas-rich blue star-forming galaxies. 

Although it is already known that the massive dark matter 


halos are dominated by massive red galaxies (e.g Zehavi et al. 
2011) and it may be inferred that those red galaxies are the 
cause of increase in clustering found by Komiya et al. (2013), 
such speculation has not been proven by any observations so 
far. This work gives the first observational evidence about the 
nature of galaxies around the AGNs as a function of BH mass 
based on statistically significant number of samples, and pro¬ 
vides an important clue to identify the feeding mechanism of 
the most massive SMBH. 

Throughout this paper, we assume a cosmology with = 
0.3, Ha = 0.7, h — 0.7 and as = 0.8. All magnitudes are given 
in the AB system. All the distances are measured in comov¬ 
ing coordinates. The correlation length is presented in unit of 
/i“^Mpc. 

2 Datasets 

2.1 AGNs 

The AGN samples were extracted from two AGN properties cat¬ 
alogs created by Shen et al. (2011) and Greene et al. (2007a). 
These catalog contains virial mass estimates of SMBHs, which 
were measured from the FWHM of emission lines and contin¬ 
uum flux density. We examined systematic difference of the 
mass estimates between the two catalogs by extracting AGNs 
which are contained in both of the catalogs. As shown in 
Figure 1, the mass estimates of Greene et al. (2007a) Mgo 7 is 
relatively lower than the estimates of Shen et al. (2011) Msn 
by 0.5 dex. We therefore calibrated Mgot with Msn using the 
following formula: 

Mgo7 = (Mgo7 - 1.06)/0.806. (1) 

We used Mqot as mass estimates of AGNs which are included 
only in the catalog of Greene et al. (2007a). For the other AGNs, 
we used Msn as their mass estimates. 

In this work the range of redshift is restricted to 0.1-1.0, 
and mass of SMBH is restricted to the range of 10® ®-10^° Mq. 
The redshift and mass distributions of AGNs used in this work 
are shown in Figures 2-A. The AGN samples from Greene 
et al. (2007a) mostly comprise lower redshift and lower mass 
samples. All the samples are divided into 12 groups shown 
in Figure 2, and analysis are performed for each group with 
enough statistics or a combined group. Sample selection is ap¬ 
plied following the criteria described in section 2.3. The number 
of analyzed AGNs is 991 from Greene et al. (2007a), 7068 from 
Shen et al. (2011), and 8059 in total. 

2.2 Galaxies 

The galaxy samples were retrieved from Virtual Observatory 
(VO) services of SDSS DR8 (Aihara et al. 2011) and UKIDSS 
DR9 LAS catalog. The UKIDSS project is defined in Lawrence 
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et al. (2007). The data retrieval and creation of the merged cat¬ 
alog were performed for each AGN. The JVO command line 
tool (Jc client ') was used to automate the data retrieval. The 
following criteria are used to retrieve the galaxy sample: For 
SDSS catalog data of which resolveStatus attrihute equals to 
257 or 258, which corresponds to select unique objects, are re¬ 
trieved. For UKIDSS LAS catalog data of which mergedClass 
attrihute equals to 1 or —3, which corresponds to select objects 
flagged as “Galaxy” or “Probable galaxy”, are retrieved. For 
both catalog data within 1 degree from the AGN coordinates 
are retrieved. 

Photometric magnitude used in this work is an aperture mag¬ 
nitude in 2 arcsec diameter for both SDSS and UKIDSS. Due to 
the difference of the point spread functions (PSF) between the 
two surveys 1.3 arcsec for SDSS in FWHM of PSF, ~0.8 arc- 
sec for UKIDSS), color measurements suffer a systematic error 
especially for objects whose apparent size exceeds the aperture 
size. In this work the comparisons of galaxy color are made 
for two samples of lower and higher BH mass, and both of the 
samples are equally affected by the error of the color estimates. 
Thus the inaccuracy in the color estimates is not so harmless to 
derive the relative difference of the distributions of galaxy color. 

The two catalogs were merged into a single one. In creat¬ 
ing the merged catalog, UKIDSS objects detected in the K band 
data were used as a reference, and a nearest neighbor was se¬ 
lected from the SDSS catalog within 2 arcsec distance for each 
UKIDSS object. The UKIDSS objects for which SDSS coun¬ 
terpart was not found were preserved in the merged catalog, 
while the SDSS objects which were not selected as a counter¬ 
part were not included in the merged catalog. Thus our galaxy 
samples are K-band selected ones. 

2.3 Data selection 

Data selection was applied to the AGN datasets to ensure their 
quality as follows: at first, coverage of the SDSS and UKIDSS 
galaxy samples were investigated by using the window func¬ 
tion for SDSS and frame metadata describing the region of ob¬ 
servation for UKIDSS. Regions outside the observed area were 
identified as a dead region. 

We found that some of the UKIDSS galaxy samples were 
spurious events caused by cosmic rays, bright stars, or unknown 
reason. The area that contains the spurious events was identi¬ 
fied by looking for a high density region in the count map of 
the samples detected in the K-band and not detected in SDSS. 
Those areas were taken as dead regions as well. 

The fractions of the dead region were calculated as a func¬ 
tion of projected distance from the AGN position with 0.2 Mpc 
bin width. If the fraction exceeds 0.2 at a projected distance less 
than 5 Mpc, the AGN dataset was removed from the samples. 

^ http://jvo.nao.ac.jp/jc_clinet/ 


Otherwise the fractions are used to correct the effective area in 
deriving the number density of galaxies. 

Next we examined the uniformity of the galaxy number den¬ 
sity around AGN. Since the number of galaxies associated with 
the AGNs is usually much smaller than the total number of 
background/foreground galaxies, a flat distribution is expected 
for the number density as a function of distance from AGN. If 
there exist a cluster/group of galaxies or stars in front of the 
AGN held, it can produce strong fluctuation in the distribution 
of number density. Inhomogeneity of the depth of observations 
can also produce discontinuous density distribution. 

To reduce a false signal produced by those cases, we calcu¬ 
lated three parameters, crmax, and Bqg, and they were used 
to filter the datasets showing non-uniformity in number den¬ 
sity. The parameter is a square sum of the deviation from 
a flat distribution, and the parameter crmax is a maximum de¬ 
viation from the average density calculated at a distance range 
from 2 to 5 Mpc. The adapted criteria are x^/(n — 1) < 3.0 and 
crmax < 5, where n is the number of data points in the galaxy 
number density distribution. Bqg was calculated as (Longair 
&Seldner 1979), 

3 — 7 Atotal — Nhg, 

^QG = -^ -, ( 2 ) 

27r po 

where Wtotai is the total number of galaxies at projected dis¬ 
tance less than 1 Mpc from AGN, Nhg is the expected number of 
background/foreground galaxies which are not associated with 
the AGN and is estimated from the number density at rp = 3.0- 
5.0 Mpc, Po is the average number density of observed galaxy at 
the AGN redshift. This parameter is usually used to estimate the 
clustering strength, but here it is used to remove datasets show¬ 
ing extraordinary large positive/negative excess density around 
AGNs. The criterion adapted is —10000 < Bqg < 10000. Only 
two samples were discarded with this selection. We checked 
that the samples for which clustering strength is smaller than 
~80 Mpc in term of correlation length don’t give rise a promi¬ 
nent fluctuation detected with the above criteria in the galaxy 
density histogram. As such, these criteria are almost free from 
selection bias caused by the clustering strength itself at the 
source. 

To maximize the signal-to-noise ratio, we removed datasets 
of shallow observations based on the parameter po, which is an 
average number density of galaxies at the AGN redshift; po is 
calculated from the luminosity function which is parametrized 
as a function of redshift z and rest-frame wavelength A as 
discussed in Komiya et al. (2013). We adapted a criteria of 
Po > Mpc“^. 

The total number of AGN datasets that have passed above 
criteria is 8059 out of the original number of 9986. About 20% 
of datasets are discarded. Among the 8059 samples, ~80% of 
the samples overlap with the Komiya et al. (2013) samples. 
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3 Analysis Method 

3.1 Correlation length 


The analysis method used for calculating the AGN-galaxy 
cross-correlation function is completely the same as that used 
in Komiya et al. (2013), which is briefly described here. 

The cross-correlation function of AGNs and galaxies ^(r) 
can be expressed as an excess in number density of galaxies 
p(r) relative to the average density po at the AGN redshift, 

«w = dr)_i. ,3) 

Po 

where r represents the distance from an AGN. We assume the 
power-law form for the cross-correlation function, 



(4) 


where ro is the correlation length and 7 is the power-law in¬ 
dex fixed to 1 . 8 , which is a canonical value measured by many 
other clustering study of galaxies and QSOs. The projected 
cross-correlation function uj{rp) is calculated by integrating 
Equation (3) as 


poo 

Lo{rp) = 2 / ^{rp,-K)dir 

Jo 



= rp 


VpJ r( 2 ) 


-0^dr 


(5) 


where tt and Vp are distance along and perpendicular to the line 
of sight, respectively, and T is the gamma function. uj(rp) can 
be derived observationally from the surface density of galaxies 
n{rp) as 


uj{rp) = 


n{rp)-nh^ 


( 6 ) 


po 

where ribg represents the surface density of back¬ 
ground/foreground galaxies which are unassociated with 
the corresponding AGN. From Equations (5) and ( 6 ), the 
surface density of galaxies around an AGN can be modeled as 


n{rp) = ( 7 ( 7 ) -po-Tp 


T ttbg, 


(7) 


where the term of the gamma function is represented hy C{'y). 
Fitting this model function to the observed surface density, we 
can obtain the best estimates of ro and ribg. As already men¬ 
tioned, 7 is fixed to 1.8 and po is calculated from the empir¬ 
ical formula of galaxy luminosity function. The details of the 
parametrized luminosity function are described in Komiya et al. 
(2013). Since the clustering signal is too weak to obtain mean¬ 
ingful parameter values for each AGN dataset, we applied the 
fitting to the average of n(rp) and po for a given AGN group. 

The uncertainty of ro is calculated as a root sum square of 
one sigma statistical error and systematic error derived from the 
uncertainty of po as discussed in Komiya et al. (2013). The 
uncertainty related with an intrinsic variance of clustering is 
also estimated from a difference between the cross correlation 


lengths calculated for two independent sub samples which are 
constructed by dividing the original samples. 

It should be noted that the cross-correlation length obtained 
by this method is not a simple average for the AGN group, but 
an average weighted with po. Thus the result is biased to the 
low-z AGN samples. 


3.2 AGN absolute bias 

The linear bias of AGNs relative to the dark matter distribution 
can be derived from the cross-correlation length between AGNs 
and galaxies as follows: Assuming linear bias, the autocorrela¬ 
tion function of AGNs ^aa is related with the cross-correlation 
function between AGNs and galaxies ^ag in the form, 

^aa = Cag/?gg, (8) 

where ^gg is the autocorrelation function of galaxies. The au¬ 
tocorrelation function of galaxies depends on galaxies color and 
luminosity, and measured by Zehavi et al. (2011) as ^gg.uiub = 
(r/3.6 and ^GG,red = (r/ 6.6 ® for 

blue and red galaxies with brightness of —20 < Mr < —19. 
The luminosity dependence is small in the range of —22 < 
Mr < —17, while for galaxies with Mr < — 22 the autocorre¬ 
lation function significantly changes to the form of ^GG.red = 
(r/10.7 /i“^Mpc)“’^'® In order to take into account the differ¬ 
ence of blue/red/brightest galaxy fraction in the samples we cal¬ 
culated the autocorrelation function of galaxies with the follow¬ 
ing formulas: 


^GG = /red^GG ,red + fbl ue Cgg ,blue ? “f/bright^GG,bright 5 (9) 


where /red, /blue, /bright IS a fractiou of red, blue, and brightest 
galaxies in the sample which are measured as described in the 
next section. The autocorrelation function of AGN is calculated 
by using Equation ( 8 ) and (9), and fitted to a power law function 
of the form of ^aa = (r/rAA 

The linear bias of AGNs b is calculated as (Koutoulidis et al. 
2013): 


( 10 ) 

where cts.agn and (T8,dm is the rms fluctuations of the AGN and 
dark matter density distribution within spheres of a comoving 
radius of 8 Mpc, respectively, cts.agn is given by 


(TS.AGN = (~^) 

> 72 ( 7 ) = 


7/2 


72 


(3-7)(4-7)(6-7)2r ’ 


( 11 ) 

( 12 ) 


and ( 78 ,DM is 


<78,DM — <78 


D{z) 

77(0)’ 


where D{z) is the linear growth factor given as. 


( 13 ) 
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Diz)- 

E{zf 


5n^E{z) 1 + y 

2 A 

■ r^m(l + z)^ + 


(14) 

(15) 


Using Equations (10), (11) and (13), bias is calculated from 
the AGN autocorrelation length as: 



/ asD{z) 

V ^(0) 


-1 


(16) 


3.3 Red galaxy fraction and luminosity distribution 

Using the UKIDSS/SDSS merged catalog obtained as described 
in section 2.2, color of each galaxy was calculated by perform¬ 
ing SED fitting. The SED fitting was performed by using EAZY 
software developed by Brammer et al. (2008). As we are inter¬ 
ested only in the galaxies associated with the AGNs, the redshift 
was fixed to the AGN redshift in the SED fitting. Photometric 
redshift was not used since the expected error can be large 
and does not improve the statistics significantly. It might also 
be possible that photo-z is not accurately determined for blue 
galaxies since the photo-z is mainly determined by the structure 
at 4000 A, so this can produce systematic bias for selecting red 
galaxies preferentially. Instead of using the photometric red¬ 
shift, we adapted filtering based on of tho SED fitting to 
select galaxies located at AGN redshifts inclusively. 

A color estimator Dopt-m is defined as 


Dopt—iR — rriopt miR, (12) 

where mopt and miR are magnitudes at wavelength range 
3,000-3,500 A and 10,000-12,000 A in the rest frame, respec¬ 
tively. In this analysis, we used only the galaxy samples which 
were detected both in UKIDSS and SDSS and for which the 
reduced of the SED fitting is less than a given limit. 

The Dopt-iR distribution for galaxies associated with the 
AGNs can be derived by the subtraction method. Dopt-in dis¬ 
tributions averaged over every AGNs in a group are calculated 
at a central and an offset region of AGN field. We defined the 
central region as the region within 1 Mpc from the AGN, and 
defined the offset region as an annulus region at a distance of 3 
to 5 Mpc. Since most of the galaxy samples at the offset region 
are located at redshifts different from those of AGNs, it can be 
regarded as background galaxies. By subtracting distribution 
at offset region from that at central region, we can obtain the 
Dopt-iR distribution for galaxies associated with the AGNs. 

The reduced x^ of tho SED fitting for galaxies at the off¬ 
set region tends to be larger than that of galaxies at the central 
region, since the fitting is performed by assuming incorrect red¬ 
shifts. We determined the limit of the reduced x^ for selecting 
galaxies, by comparing the distributions for central and offset 
galaxies. We set the limit to 4.5 for z < 0.6, and 2.0 for z > 0.6. 

In Figure 5, we show the distribution of Dopt-in derived 


from all the AGN samples. The bimodality of galaxy color dis¬ 
tribution is clearly seen. In the same figure, the color distri¬ 
butions for blue and red galaxies are deconvolved by assuming 
normal distribution. The peak values of the color distribution 
for blue and red galaxies are 3.1 and 4.2, respectively. We de¬ 
rive the red galaxy fraction for each AGN group by fitting the 
color distribution to a model function given with a sum of two 
normal distributions. 

The luminosity distribution can be derived with the same 
method as that used for the color distribution described above. 
The absolute magnitude was calculated by miR — DM{z), 
where DM (z) is the distance modulus at redshift z. In Figure 6, 
we show the color-magnitude diagram obtained for all the AGN 
samples by the subtraction method. The red sequence and blue 
cloud are clearly seen. Red sequence galaxies are concentrated 
in a upper right region (M < —20 and Dopt-in > 3.5), while 
blue cloud galaxies are widely spread in the dimmer part. 


4 Results 

4.1 Cross-correlation 

All the AGN samples are divided into 12 redshift-mass groups 
as shown in Figure 2. Redshift is divided into three ranges 
0.1-0.3, 0.3-0.6 and 0.6-1.0. We designate them as zl, z2, 
and z3, respectively. BH mass is divided into four ranges 
log(MBH/Mo) = 6.5-7.5, 7.5-8.2, 8.2-9.0, 9.0-10.0, and 
they are designated as M65, M75, M82, and M90, respectively. 
Hereafter we call, for an example, a group of redshift range of 
0.1-0.3 and mass range of 6.5-7.5 as M65-zl. 

The results obtained by fitting the model function of 
Equation (7) to the observed average number densities are sum¬ 
marized in the first part of Table 1. In Figure 7, the cross¬ 
correlation length ro obtained for each redshift-mass group is 
plotted as a function of Mbh- The result shows that there is a 
tendency that ro increases with Mbh, while there is no clear 
redshift dependence. 

As no clear redshift dependence is seen, the AGN samples 
are combined for the same mass and all the redshift groups to 
reduce the uncertainty resulting from intrinsic variance of clus¬ 
tering around each AGN, then ro is calculated for the combined 
mass groups. The projected number densities for these samples 
and fitted model functions are shown in Figure 8. The horizon¬ 
tal dashed lines in these figures correspond to nog obtained from 
the fitting. The projected cross-correlation functions calculated 
from equation (6) are shown in Figure 9, and are well expressed 
by a power law function. 

The cross-correlation lengths derived from the fitting are 
plotted in Figure 10. The error bars include statistical error, 
systematic error, and uncertainty related with the intrinsic vari¬ 
ance among AGNs. The uncertainty due to the intrinsic vari¬ 
ance was estimated from a difference between the cross corre- 
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lation lengths calculated for two independent sub samples. The 
sub samples were constructed by dividing the original AGN 
samples into two groups in random. We generated 21 sets of 
two half-sized samples, calculated absolute differences of cross 
correlation lengths between the two sub samples ( 5 ro,i /2 for 
each set, and selected median of 5ro 1/2 as an estimate of the 
standard deviation for the half-sized sub sample. The corre¬ 
sponding error on ro of the original sample was calculated as 
Org = which is summarized in Table 2. 

The difference of ro between M75 and M90 groups is about 
2.6 sigma excluding the systematic errors. For comparison the 
results obtained in the previous work (Komiya et al. 2013) are 
also plotted in the figure. They are consistent with the results of 
this work. 

The AGN linear bias is estimated from the cross-correlation 
length assuming that the autocorrelation function of the galax¬ 
ies is expressed by a linear combination of autocorrelation func¬ 
tions of three types of galaxies. 

Zehavi et al. (2011) measured an autocorrelation function 
of galaxies up to redshift of 0.25, and derived its color and 
luminosity dependence. Their results show that the auto¬ 
correlation function significantly different among blue, red, 
and the brightest galaxies. Referring to their results, we 
used autocorrelation functions of the form CGG,biue(r) = 
(r/3.6h“^Mpc)“^'^, CGG.red(r) = (r/6.6h“^Mpc)“^ ®, and 
^GG.bright (r) = (r/10.7h,“’^Mpc)“^ ® for blue and red galaxies 
with absolute brightness in r band Mr > —22, and the bright¬ 
est galaxies with Mr < —22, respectively. The parameter val¬ 
ues adapted for blue and red galaxies are the ones obtained for 
galaxies with —20 < Mr < —19. Although the absolute magni¬ 
tude of our galaxy samples ranges over — 23 < Mr < —17, we 
ignore the luminosity dependence in the range of — 22 < Mr < 
— 17 as their measured dependence is small. For galaxies with 
brighter than Mr = — 22 significant increase in autocorrelation 
function was observed, so we used a separate function for those 
galaxies. We also assume that the autocorrelation functions are 
unchanged at redshifts up to 1 . 0 . 

To derive weighting factors among these three components, 
fractions of blue and red galaxies with Mr > — 22 and of galax¬ 
ies with Mr < —22 are measured by the subtraction method de¬ 
scribed in section 3.3. The obtained fractions are summarized 
in Table 2. 

Using the autocorrelation function of galaxies derived from a 
linear combination of those for three galaxy types, we obtained 
b= 1.52^0.54, 1.471q 38, 2.38tQ;49, for mass group 

M65, M75, M82 and M90, respectively. The detailed method 
for deriving the AGN bias is described in section 3.2. These 
results are summarized in Table 2 and plotted in Figure 11. The 
uncertainty of b includes statistical error, systematic error cor¬ 
responding to the uncertainty of po estimation, and the intrinsic 
variance of the clustering around each AGN. 


As is inferred from the mass-redshift distribution shown in 
Figure 2, even for the same redshift group the higher mass 
groups tends to be biased to higher redshift. To reduce the red¬ 
shift bias, we constructed redshift matched samples for each 
redshift group. The redshift matched samples were constructed 
by selecting the same number of AGN samples for each redshift 
bin with 0.02 width. 

The number of samples selected in this way is reduced sig¬ 
nificantly from the original samples, so the intrinsic fluctuation 
of clustering strength of each sample can introduce large scat¬ 
ter to the cross-correlation length. To obtain the most typical 
sample in term of clustering, we randomly constructed 21 sets 
of redshift matched sample for each redshift group, calculated 
cross-correlation length for each sample set, and selected a sam¬ 
ple set for which cross-correlation length was median among 
the sample sets. The result is shown in Figure 12 and summa¬ 
rized in the middle part of Table 1. The tendency of increase in 
ro as a function Mbh is seen in the redshift matched samples as 
well. 

To see the redshift dependence of ro for the same Mbh, we 
constructed mass matched samples in the same way as used to 
construct redshift matched samples with 0.2 dex of bin width of 
mass histogram. The results are shown in Figure 13 and sum¬ 
marized in the last part of Table 1. No clear redshift dependence 
is seen for any mass group. 

These results are consistent with our previous result in 
Komiya et al. (2013), although we have used the most recent 
UKIDSS catalog and slightly different data selection criteria 
and mass correction. 

4.2 Properties of Galaxies 

To investigate the properties of galaxies at environment of 
AGNs with the most massive SMBHs, where increase of clus¬ 
tering is found, we calculated the color parameter Dopt-m for 
all the detected galaxies and derived its distributions for galax¬ 
ies associated with the AGNs by using the subtraction method 
described in section 3.3. Dopt-ia parameter corresponds to 
the color defined as difference in brightness between the wave¬ 
length ranges of 3,000-3,500 A and 10 , 000 - 12,000 A in the rest 
frame. In Figure 14 we show Dopt-m parameter distributions 
for each redshift and mass sample. The observed distributions 
are fitted with a two component model (red and blue galaxy 
components) assuming the normal distribution for both compo¬ 
nents. The mean and standard deviation of the distribution of 
each component are fixed to the predetermined values, and only 
the mixing ratio and the normalization constant are taken as free 
parameters. The predetermined values of the mean and standard 
deviation are obtained by fitting the model function to the ob¬ 
served distribution for all the AGN samples with taking all the 
parameters free. The parameter values obtained from this fitting 
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are summarized in Table 3. 

The fractions of red galaxies are plotted in the left panel of 
Figure 15 as a function of Mbh- The normalized excess densi¬ 
ties Lu', which are defined as lu' = (rion —noff)/po, are plotted in 
the right panel of the figure, non and noff are the surface num¬ 
ber density of galaxies at projected distance of 0-lMpc (on re¬ 
gion) and 3-5Mpc (off region) from AGN, respectively. The oj' 
can be considered as an approximate estimate of the projected 
cross-correlation function. The error bars represent one sigma 
Poisson statistical error for both figures. Since red galaxies are 
typically brighter than blue galaxies as seen in Figure 6, the ob¬ 
served red galaxy fraction tends to be higher for the higher red- 
shift samples. Thus, the comparison of the red galaxy fraction 
is meaningful only for the same redshift samples. 

At redshift zl, both of the red galaxy fraction and the nor¬ 
malized excess density are almost unchanged among the the 
mass groups of M65, M75 and M82. At redshift z2, significant 
increase of the normalized excess density is seen for the red 
galaxies of the mass group M90. At redshift z3, blue galaxies 
are hardly detected and red fraction is ~100% for mass groups 
M82 and M90. 

The data samples used in Figures 14 and 15 have a slight 
difference in the redshift distribution among the respective mass 
groups even for the same redshift groups; the lower mass group 
has relatively smaller redshift than the higher mass group. The 
red fraction of the higher mass group, therefore, can be biased 
to higher value. To reduce the effect of redshift bias, the same 
analysis was performed also for the redshift matched samples. 

The distributions of Dopt-ia for the redshift matched sam¬ 
ples are shown in Figure 16. Since the statistics of lower mass 
groups are poor, the two lower mass groups are combined. The 
obtained red galaxy fractions and normalized excess density are 
plotted as a function of Mbh in Figure 17. These results also 
show that red galaxy becomes the dominant component in the 
highest mass group, M90. 

To see the relative difference of luminosity function of galax¬ 
ies around SMBHs with lower and higher mass, absolute magni¬ 
tude distributions are compared in Figure 18. They are obtained 
by using the subtraction method described in section 3.3. In this 
comparison, redshift matched samples are used so that both of 
the samples have the same sensitivity on the detection of galax¬ 
ies. To check the equivalence among the samples in term of the 
sensitivity of observations, the distributions of po are compared 
in Figure 19. We found that they are consistent with each other. 

The left panel of Figure 18 shows the comparison between 
the AGN groups of M65-I-M75 and M82 for redshift range of 
zl, and the ratio (high mass/low mass) of the magnitude distri¬ 
butions is shown in the bottom of the panel. Since the lower 
mass samples has smaller excess density than the higher mass 
sample, we combined two lower mass samples to increase the 
statistics. The red histogram is the distribution of absolute mag¬ 


nitude for higher mass sample, and the blue histogram is for 
lower mass sample. The absolute magnitude is estimated at 
wavelength range of 1-1. 2A in the rest frame by fitting model 
SED to the observed one using the EAZY code as described 
in section 3.3. The ratio is peaked around —20.25 mag, and it 
shows a steep decline at dimmer side. The significance of the 
depletion of galaxies in the M82-zl group against (M65-I-M75)- 
zl is ~2.5 sigma at magnitude range from —19.5 to —18.0 mag. 
The dashed line on the ratio plot is a power law function fitted 
to the observations at brighter side below —20.0 mag. The best 
fit value of the power law index is —0.096±0.075 for the ratio 
between M82-zl and (M65-i-M75)-zl samples. 

The same comparison between the mass groups M75-I-M82 
and M90 at redshift z2 and z3 is shown in the right panel of 
Figure 18. The observed ratio is fitted with the power law func¬ 
tion and the power law index is estimated to be -1-0.046 ± 0.093, 
which is consistent with a constant over the luminosity range. 
The ratio of M82 and M65-I-M75 at redshift zl and that of M90 
and M75-I-M82 at redshift z2 and z3 are compared in Figure 20. 

5 Discussion 

We have successfully reconfirmed the increase of cross¬ 
correlation length above Mbh = 10®'^ Mq (Figure 10), that 
was already reported in the previous paper by Komiya et al. 
(2013), by adapting the improved data selection and using the 
most recent UKIDSS DR9 dataset. We also confirmed that the 
trend is also seen in the redshift matched samples (Figure 12), 
and the cross-correlation length does not depend on the redshift 
(Figure 13). Although the AGN luminosity dependence was not 
investigated in this paper, the previous work reported no signif¬ 
icant dependence on the luminosity (Komiya et al. 2013). 

The derived AGN bias b (Figure 11) increase with Mbh 
above 10®'^ Mq. The AGN bias for M75 mass group is b = 
1.471q 3 ® at average redshift of 0.47, which is consistent with 
those of SDSS QSOs by Ross et al. (2009), 2dF QSOs by Groom 
et al. (2005), X-ray selected AGNs by Hickox et al. (2009), and 
optically selected AGNs by Krumpe et al. (2012). 

The AGN bias of mass group M90 is estimated to be b = 
S.OS^Q g®, which is larger than that of radio AGNs of Hickox et 
al. (2009) by one sigma and indicates that those AGNs reside in 
haloes with mass larger than Mq. 

Shen et al. (2009) also studied virial mass dependence of 
clustering strength using SDSS QSOs at redshift range from 0.4 
to 2.5. They found that the difference in the clustering strength 
for the 10% most massive QSOs and the remaining 90% is sig¬ 
nificant at the ~ 2(7 level. They also studied radio activity de¬ 
pendence of clustering by comparing between the radio-loud 
and radio-quiet QSOs samples, and found that the radio-loud 
QSOs are more strongly clustered than the radio-quiet QSOs at 
the ~ 2.5 ct level. Although these radio-loud QSOs tend to have 
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systematically larger BH masses than those radio-quiet QSOs 
by 0.12 dex, the difference in their clustering remains when the 
comparison is made for mass matched samples. Thus they ar¬ 
gue that more massive host halos, and denser environments may 
be related to the triggering of radio activity. 

Mandelbaum et al. (2009) also found that the radio AGNs 
are hosted by dark matter halo with ~1.6x 10^®h“^ Mq. There 
are observational evidences that radio-loud AGNs are associ¬ 
ated with the massive galaxies (Best et al. 2005) and have 
SMBHs with masses typically larger than 10® Mq (Laor 2000). 
Ishibashi et al. (2014) also claimed that radio galaxies are pref¬ 
erentially associated with the more massive black holes. 

According to the theoretical predictions by Fanidakis et al. 
(2011) for radio loudness and black holes mass (see bottom 
panel of Figure 14 in their paper), AGNs in the hot halo mode 
are distributed at a region of higher radio-loudness and their 
black hole mass extends from 10® Mq up to lO’^® Mq, while 
the AGNs in the cold accretion mode are at lower part in radio¬ 
loudness and their black hole mass is limited below ~ 10® '® Mq . 
It is, therefore, expected that AGNs in the hot halo mode, i.e. 
in dark matter haloes with higher mass, can preferentially be 
selected by higher radio-loudness or higher black home mass. 
Thus the increasing trend of AGN bias at higher black hole 
masses found in this work are consistent with those results ob¬ 
tained by the other authors, and they all are also consistent with 
the picture drawn from Fanidakis et al. (2011). 

The main purpose of this paper is to investigate what type 
of galaxy contributes to the increase of galaxy density around 
the most massive SMBH. In this work we considered two kind 
of galaxy types which are classified based on the distribution in 
a color-magnitude diagram. One is a blue cloud galaxy (blue 
galaxy) which occupies a bluer and dimmer side of the distribu¬ 
tion, and the other is a red sequence galaxy (red galaxy) which 
has a redder and narrower distribution in the color-magnitude 
diagram (Figure 6). 

It should be noticed that detection efficiency for blue galaxy 
decreases more rapidly with increasing redshift than for red 
galaxy due to the difference in their brightness. As is shown in 
the bottom panels of Figure 14, at redshift range of z3 (z =0.6- 
1 .0) the observed fraction of blue galaxies is very low, this is 
largely because most of the blue galaxies are below the detec¬ 
tion limit. At redshift range of z2 {z =0.3-0.6) the fraction of 
the detected blue galaxy constitutes about 20% of all the de¬ 
tected galaxies, while at redshift range of zl {z =0.1-0.3) the 
fraction increase to around 40%. As shown here, the observed 
red/blue fraction has strong redshift dependence and should not 
be compared between different redshift groups. 

The results of the two component analysis on the Dopt-m 
distribution (Figures 14,15) indicate that the increase of cluster¬ 
ing strength found in the cross-correlation analysis for the most 
massive mass group M90 relative to the lower mass groups are 


mainly due to the contribution of red galaxies. This can be jus¬ 
tified from the observational evidence that increase of the nor¬ 
malized excess density at the transition from M82 to M90 for 
redshift z2 is significant for red galaxies, while it is almost un¬ 
changed or rather decreasing for blue galaxies. It is known that 
early type galaxies tend to be found at a high density region 
such as a cluster/group core (Dressier 1980; Postman & Geller 
1984; Balogh et al. 2004). It is, therefore, naturally expected 
that the high density around the high mass SMBH is coupled 
with early type galaxies. 

The other indication from this analysis is that the increase 
of the clustering strength found in the cross-correlation analy¬ 
sis for M82 relative to the M75 is due to the contribution from 
both of the red and blue galaxies. The red fraction is almost 
unchanged between M82 and M75, and slight increase in the 
normalized excess density from M75 to M82 is seen in both of 
the red and blue components for redshift groups zl and z2. 

The same analysis is performed also for the redshift matched 
samples (Figure 16 and 17). The trend of the increase in the 
red fraction and normalized excess density for red galaxies of 
the most massive group (M90) is also seen in this result. The 
comparison between M82 and M65-I-M75 at redshift zl shows 
that increase in the normalized excess density are seen in both 
the red and blue galaxies. 

Two physical mechanisms could be relevant to this evolu¬ 
tion. One is related with the increase of brightness of blue dim 
galaxies induced by starburst triggered in a secular and merger 
mode, which can result in the increase of the detectable galaxy 
density. Another one is related with the transition from blue to 
red galaxies, which can be the result of a feedback or gas re¬ 
moval mechanism shutting off the star burst activity. The latter 
process may follows the former process. In case where there 
is an enough time delay between the both processes, we could 
observe the increase of density both for red and blue galaxies as 
is the case of M82 group. If the time delay is short, only the in¬ 
crease of red galaxy density will be observed, and this could be 
the case of M90 group. The evolution of properties of galaxies 
around AGNs observed in this work could be explained in this 
manner. 

According to the theoretical works by Fanidakis et al. 
(2013), the hot halo mode can be the most effective accretion 
mechanism in massive haloes with Mh > Mq. In such 

a massive halo most of the galaxies are expected to be in red 
sequence since AGN feedback shut off the gas cooling and star 
formation in their model. In galaxies turned into red sequence, 
i.e. early type galaxies, gas accretion from the evolving stars 
to the center of SMBH could also occur as claimed by Ciotti & 
Ostriker (2001). Considering these theoretical predictions, the 
dominance of red galaxies around AGNs with higher BH mass 
obtained in this work is indicative of predominance of hot halo 
mode and/or gas accretion from evolving stars for the growth of 



10 


Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0 


the most massive SMBHs. 

The Eddington ratios of our AGN samples in M90-zl group 
are distributed in the range A = 0.003 - 0.1, and peaked around 
log(A) ~ —1.8. The lower bound of the distribution is limited 
by the detection limit, so the intrinsic peak can be expected to 
be much lower than log(A) = —2. The bolometric luminosity of 
the same samples is above lO'^® erg/s. According to the model 
of Fanidakis et al. (2013b) for those bolometric luminosity, two 
branches exist in the Mhaio-Tboi diagram (see Figure 1 in their 
paper); one is the branch of starburst mode for AGNs in dark 
matter halo of 10^^ Mq, the other is the branch of hot halo 
mode for AGNs in dark matter halo of 10^® Mq. Considering 
that the Eddington ratios of our samples are lower than those 
expected for the starburst mode (~0.1), the sample of M90-zl 
are consistent with the AGNs in the branch of hot halo mode. 

We also investigated relative difference of luminosity func¬ 
tion between the low and high mass groups. The redshift 
matched samples are used for making them directly comparable 
with each other. The equivalence among the samples in term of 
the sensitivity of observations was checked by comparing the 
distributions of po (Figure 19), and we found that they are con¬ 
sistent with each other. Thus the distributions of completeness 
fraction are the same for the redshift matched samples, and the 
completeness fraction can be canceled by calculating the ratio. 

For the redshift range zl, comparison between the M82 mass 
sample and the lower mass sample, which is a combined sample 
of M65-I-M75, is made in the left panel of Figure 18. The figure 
indicates that an increase of galaxy density for the higher mass 
sample would be caused by galaxies brighter than Mir = —19.5 
mag. At a dimmer side of Mir > —19.5 mag, on the other hand, 
the luminosity function of the higher mass sample is smaller 
than that of the smaller mass group. The significance of fhe 
deficit of galaxies for the higher mass sample relative to the 
smaller mass sample is 2.5 sigma at absolute brightness range 
from —19.5 to —18.0 mag. 

Considering that blue galaxies dominate the dimmer side of 
the luminosity function, the difference of the luminosity func¬ 
tions could be the result of brightening of the dim blue galaxies 
in the transition from lower mass group (M65+M75) to higher 
mass group (M82). Besides, to make the red/blue galaxy ratio 
unchanged some but not all of the blue galaxies are need to be 
transformed to red galaxies. 

The plot of the ratio between luminosity function of M82 
and M65+M75 shows that the increase fraction is larger at the 
dimmer side and peaked around Mir = —20.5 mag, and it is 
approximately fitted with a power law function with the power 
index of —0.096±0.075. This may indicate that the dimmer 
galaxies have higher probability to increase their brightness by 
interaction and/or major/minor merger. 

For the redshift ranges z2 and z3, a comparison between 
the M90 and M75-I-M82 groups is made in the right panel of 


Figure 18, which shows the ratio of the luminosity functions is 
almost constant at Mir < —20 mag. If the brightening of blue 
galaxies occurs at a same rate between the two groups, the dif¬ 
ference is only the rate of blue to red transformation, and the 
transformation occurs without much affecting the brightness in 
the 1-1.2 A band, the constant ratio of the luminosity functions 
is naturally expected. 

6 Conclusions 

In this paper, using the updated UKIDSS catalog we have suc¬ 
cessfully reconfirmed the previous findings of Komiya et al. 
(2013) that the clustering of galaxies around AGNs with the 
most massive SMBH is larger than those with less massive 
SMBH. The AGN bias was derived for each BH mass group. 
The obtained AGN bias are fo= 1.521 q®4, 1 .47lg'3g, 2.38lQ4g, 
S.OSIq gj for mass group M65, M75, M82 and M90, respec¬ 
tively. We further investigated what type of galaxies are associ¬ 
ated in the environment of the most massive SMBH to reveal the 
nature of evolution mechanism of SMBHs and galaxies. As a 
result, it is found that red galaxies are dominated around AGNs 
with the most massive SMBH. The red galaxy fraction increases 
from 0.73±0.04 for M82 mass group to 0.98±0.07 for M90 at 
redshift z2. This is the first observational result that revealed the 
nature of galaxies at such environments from more than thou¬ 
sands of AGN samples. We also compared, for the first time, 
luminosity functions of galaxies at environments of AGNs with 
lower and higher mass SMBH, and found that there is an indica¬ 
tion that shows brightening of dim and blue galaxies along with 
the evolution of SMBHs. Summarizing the results obtained in 
this work, we can deduce the following scenario on the evo¬ 
lution of SMBH and its environment galaxies as a function of 
Mbh. 

Below the critical mass around 10®'^ Mq, the environment 
of SMBH does not depend on Mbh and is almost equivalent 
to that of a quiescent galaxy. At that environment, the fraction 
of blue galaxy is more than 40%. Secular evolution can be the 
main driver for the evolution of these SMBHs. Above the criti¬ 
cal mass, dim and blue galaxies increase its brightness, which is 
presumably caused via starburst. Some of the brightened blue 
galaxies are transformed to red galaxies after the AGN feedback 
or some other mechanisms which shut off star formation in the 
galaxies operates. 

In the environment of AGNs of M82 mass group, increases 
in densities for both blue and red galaxies are observed, while 
in the environment of AGNs of M90 mass group, the increase 
is seen only for the red galaxies. This difference may be due 
to the difference of time lag between the starburst and some 
mechanism to shut off the starburst activity. 

Considering that there is a correlation between the transition 
in properties of environmental galaxies and mass evolution of 
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SMBHs above Mbh = 10®'^ M©, they are governed by a sort of 
environmental effect rather than internal secular evolution. This 
transition makes galaxies observable by pushing up the bright¬ 
ness above the detection limit, and results in the increase of 
observed galaxy number density and the clustering length. 

The hot halo mode as well as gas accretion from evolving 
stars may be one of the most dominant mechanism in the growth 
of the most massive SMBHs, as inferred from properties of en¬ 
vironment galaxies. 
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Table 1. Parameters derived from the cross-correlation analysis for each BH mass and redshift range 


BH mass 

redshift 

T)^ 

'^AGN 

(logM) 


rd 

“^stat 

(ribg)^ 

(po)® 

logM 





h~^Mpc h~ 

+ Mpc 

Mpc-2 

10-3Mpc-3 

All samples 

6.5-7.5 

O.l-l.O 

524 

7.2 

0.24 

'^•-‘--0.6 

0.34 

38.85±0.05 

5.08±0.8 


0.1-0.3 

429 

7.1 

0.20 

'^•^-0.6 

0.36 

45.17±0.06 

5.61±0.9 

7.5-8.2 

O.l-l.O 

1880 

7.9 

0.47 

c; Q+U.S 
^•^-0.6 

0.23 

15.24±0.02 

2.71±0.5 


0.1-0.3 

520 

7.8 

0.22 

^•'^-0.6 

0.32 

37.60±0.05 

5.32±0.9 


0.3-0.6 

804 

7.9 

0.45 

4 q+0.8 

0.35 

9.41 ±0.02 

2.48±0.5 


0.6-1.0 

556 

8.0 

0.75 

K 4 + 1.6 

0.^1 2 

0.85 

2.78±0.01 

0.61±0.2 

8.2-9.0 

O.l-l.O 

4485 

8.6 

0.65 

7 7+1.^ 

'•'- 0.8 

0.16 

6.89±0.01 

1.42±0.3 


0.1-0.3 

362 

8.5 

0.23 

7 7+0.9 

'•'-0.7 

0.32 

33.10±0.06 

5.18±0.8 


0.3-0.6 

1380 

8.6 

0.46 


0.21 

8.59±0.02 

2.29±0.5 


0.6-1.0 

2743 

8.6 

0.80 


0.33 

2.58±0.01 

0.49±0.2 

9.0-10.0 

O.l-l.O 

1170 

9.2 

0.72 

g.otli 

0.34 

4.79±0.01 

1.01±0.3 


0.3-0.6 

289 

9.2 

0.48 


0.42 

8.19±0.04 

2.17±0.5 


0.6-1.0 

845 

9.2 

0.83 

Q y + 2.8 

y. t i 7 

0.62 

2.50±0.01 

0.44±0.2 

Redshift matched samples 

6.5-7.5 

0.1-0.3 

298 

7.1 

0.22 

'^•'^-0.6 

0.44 

35.04±0.07 

5.29±0.8 

7.5-8.2 


298 

7.8 

0.22 

6 3+'^-® 

^•'^-o.e 

0.41 

36.23±0.07 

5.35±0.9 

8.2-9.0 


298 

8.5 

0.22 

o 1 +0.9 

0.34 

35.96±0.07 

5.33±0.9 

7.5-8.2 

0.3-0.6 

289 

7.9 

0.48 

^•'^-0.8 

0.65 

8.47±0.03 

2.18±0.5 

8.2-9.0 


289 

8.6 

0.48 

6-9lo,9 

0.49 

8.25±0.03 

2.17±0.5 

9.0-10.0 


289 

9.2 

0.48 

8-8ll:^ 

0.42 

8.19±0.04 

2.17±0.5 

7.5-8.2 

0.6-1.0 

483 

8.0 

0.76 

5 0 + 1.7 

0.97 

2.73±0.01 

0.58±0.2 

8.2-9.0 


483 

8.6 

0.76 

o 1 + 2.2 
o.l_1.4 

0.72 

2.67±0.02 

0.58±0.2 

9.0-10.0 


483 

9.2 

0.77 

Q 2 + 2.5 

0.68 

2.76±0.02 

0.58±0.2 

Mass matched samples 

7.5-8.2 

0.1-0.3 

371 

7.9 

0.22 

6 

'^•-'--0.6 

0.38 

36.91±0.06 

5.30±0.9 


0.3-0.6 

371 

7.9 

0.45 

4 7+0.8 

^•'-0.7 

0.54 

9.53±0.03 

2.45±0.5 


0.6-1.0 

371 

7.9 

0.74 

4 4+1-0 

1.17 

2.85±0.02 

0.64±0.2 

8.2-9.0 

0.1-0.3 

362 

8.5 

0.23 

1 + 
0 0 

0.32 

33.10±0.06 

5.18±0.8 


0.3-0.6 

362 

8.5 

0.46 

6-8loi 

0.42 

8.54±0.03 

2.32±0.5 


0.6-1.0 

362 

8.5 

0.80 

7 7-I-2.3 

'-1.5 

0.99 

2.61±0.02 

0.50±0.2 

9.0-10.0 

0.3-0.6 

289 

9.2 

0.48 


0.42 

8.19±0.04 

2.17±0.5 


0.6-1.0 

289 

9.2 

0.84 

9.4+3° 

1.18 

2.49±0.02 

0.40±0.2 


“number of sample AGNs. 

^average of logarithm of BH mass. 

“average redshift. 

‘^correlation length, the error contains systematic eiTor due to uncertainty of po and Icr statistical error. 
“ statistical eiTor of tq . 

^ average of projected number density of background galaxies. 

® average of the averaged number density of galaxies at the AGN redshift. 


Table 2. Parameters derived from the bias estimate 


mass 

redshift 

fa 

^ blue 

fh 

red 

fC 

bright 

^ro 

^AA 

+A 

bs 

log Mq 






1 

n 

1 

n 



6.5-7.5 

0 . 1 - 

- 1.0 

0.321 

0.669 

0.010 

1.16 

6.33tli 

1.69 

C 9 + O .66 
J--OZ-0.54 

7.5-8.2 

0 . 1 - 

- 1.0 

0.296 

0.650 

0.054 

0.65 

5.29+2+ 

1.69 

1 47+0.49 
-'■^'-0.38 

8.2-9.0 

0 . 1 - 

- 1.0 

0.230 

0.581 

0.189 

0.35 

8.45t3;2 

1.69 

9 ‘^S+^- ' ^ 

^.00_Q 4g 

9.0-10.0 

0 . 1 - 

- 1.0 

0.112 

0.668 

0.220 

0.82 

11.09l®+ 

1.68 

3.08tJ;2? 


“fraction of blue galaxies with brightness Mr > —22 
^fraction of red galaxies with brightness Mr > —22 
“fraction of galaxies with brightness Mr < —22 

‘^uncertainty of the cross-correlation length tq due to the intrinsic variance of clustering around each AGN. 
“autocorrelation length of AGNs. 

^power law index of AGN autocorrelation function. 

®AGN bias. 
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Table 3. Parameters derived from the two component analysis on -Dopt-m distributions 


redshift 

mass 

logM© 

/^red 

^red 

/"blue 

^blue 


/ / f 

^red 

lO^Mpc 

^blue 

lO^Mpc 

All samples 

O.l-l.O 

6.5-10.0 

4.17±0.05 

0.41±0.03 3.06±0.23 

0.50±0.13 

0.77±0.09 



0.1-0.3 

6.5-7.5 

4.17 

0.41 

3.06 

0.50 

0.66±0.06 

0.64±0.09 

0.32±0.06 


7.5-8.2 

4.17 

0.41 

3.06 

0.50 

0.66±0.05 

0.60±0.07 

0.31 ±0.05 


8.2-9.0 

4.17 

0.41 

3.06 

0.50 

0.66±0.05 

0.73±0.08 

0.37±0.06 

0.3-0.6 

7.5-8.2 

4.17 

0.41 

3.06 

0.50 

0.84±0.10 

0.43±0.08 

0.08±0.06 


8.2-9.0 

4.17 

0.41 

3.06 

0.50 

0.73±0.04 

0.67±0.06 

0.25±0.05 


9.0-10.0 

4.17 

0.41 

3.06 

0.50 

0.98±0.07 

1.37±0.15 

0.22±0.10 

0.6-1.0 

7.5-8.2 

4.17 

0.41 

3.06 

0.50 

0.81±0.30 

0.40±0.17 

0.09±0.17 


8.2-9.0 

4.17 

0.41 

3.06 

0.50 

1.00±0.08 

1.28±0.09 

0.00±0.10 


9.0-10.0 

4.17 

0.41 

3.06 

0.50 

1.00±0.14 

1.39±0.17 

0.00±0.19 

Redshift matched samples 

0.1-0.3 

6.5-8.2 

4.17 

0.41 

3.06 

0.50 

0.64±0.04 

0.73±0.07 

0.41±0.05 


8.2-10.0 

4.17 

0.41 

3.06 

0.50 

0.66±0.05 

0.92±0.10 

0.48±0.08 

0.3-0.6 

7.5-9.0 

4.17 

0.41 

3.06 

0.50 

0.59±0.09 

0.41±0.10 

0.28±0.07 


9.0-10.0 

4.17 

0.41 

3.06 

0.50 

0.98±0.07 

1.37±0.15 

0.02±0.10 

0.6-1.0 

7.5-9.0 

4.17 

0.41 

3.06 

0.50 

0.86±0.16 

0.65±0.14 

0.11±0.14 


9.0-10.0 

4.17 

0.41 

3.06 

0.50 

1.00±0.13 

1.50±0.20 

0.00±0.19 


^mean of the .Dopt—IR distribution for red component. 

^standard deviation of the i!)opt —IR distribution for red component, 
^mean of the i^opt —IR distribution for blue component. 

^standard deviation of the i^opt —IR distribution for blue component, 
^fraction of red component. 

^normalized excess density for red component. 

^noiTTialized excess density for blue component. 



log(MBH/M,un) (Shen+2011) 

Fig. 1. Comparison of black hole masses derived by Shen et al. (2011) (horizontal axis) and Greene et al. (2007a) (vertical axis). The solid red line represents 
a liner function fitted to the data points. 
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2. Distribution of redshift and BH mass of AGNs. AGNs from Greene et al. (2007a) and Shen et ai. (2011) are plotted in red and blue colors, respectively, 
mass correction is applied to AGNs of Greene et al. (2007a) to compensate the systematic difference between the two catalogs. 



log(MBH/Msoi) 


Fig. 3. Distribution of BH mass for the same sample as shown in Figure 2. The red histogram with smaller masses is for Greene et al. (2007a), the blue one is 
for Shen et al. (2011). The sum of them is shown with the black histogram. 
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Fig. 4. Distribution of redshift for the same sampie as shown in Figure 2. The red histogram with smaiier redshift is for Greene et ai. (2007a), the blue one is 
for Shen et al. (2011). The sum of them is shown with the black histogram. 


-3 



Fig. 5. TTopt-iR parameter distribution for ail the AGN samples is shown with filled black circles with statistical error bars. The distribution is fitted with a 
two Gaussian model (black solid line). Each component represented with a single Gaussian corresponds to blue cloud (peak with smaller Z>opt-iR) and red 
sequence galaxies (peak with larger Uopt-m). 
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Absolute Mag (l-1.2unn) 

Fig. 6. Color magnitude distribution for ail the AGN samples. Densities are represented with gray scales and contours of equi-density. Red sequence galaxies 
are concentrated in a upper right region (M < -20 mag and £)opt-iR > 3.5), while blue cloud galaxies are widely spread in the dimmer part. 


U 

Q. 



6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 


log(MBH/Msoi) 


Fig. 7. Cross-correlation lengths derived for each redshift and mass group are plotted as a function BH mass. The data point of the same redshift groups are 
connected with a line. Solid circles are for zl redshift group, solid triangles and boxes are for z 2 and z3 redshift groups, respectively. 
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Projected distance from AGN [Mpc] Projected distance from AGN [Mpc] 




Projected distance from AGN [Mpc] Projected distance from AGN [Mpc] 

Fig. 8. Surface densities of galaxies plotted as a function of projected distance from AGN for mass group M65, M75, M82 and M90. 



Tp (Mpc) 



Tp (Mpc) 



Tp (Mpc) 


Fig. 9. Projected cross-correlation functions for mass group M65, M75, M82 and M90. 
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14 


12 


10 

Q. 

8 

I 

O 

^ 6 


4 


2 

6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 

log(MBH/Msoi) 

Fig. 10. Cross-correlation lengths derived for each mass groups are plotted as a function of BH mass. The redshift range of each mass groups are 0.1~0.3 
for M65 mass group, 0.1 .0 for M75 and M82 mass groups, and 0.3~1.0 for M90 mass group. Results of this work are shown with filled circles and previous 

work by Komiya et al. (2013) are shown with open circles. The error bars of this work include statistical error, systematic error, and uncertainty related with the 
intrinsic variance among AGNs. 


z = 0.1 ~ 1.0 
—This work 




0.0 0.5 1.0 1.5 2.0 2.5 


Redshift 

Fig. 11. AGN bias plotted as a function of redshift. The results of this work are plotted with filled circles for mass groups of M65, M75, M82 and M90. The 
resulfs by Hickox et al. for radio. X-ray and IR selected AGNs are plotted with open squares at b = 2.03, 1.40, and 0.95, respectively. The results by Krmpe et 
al. for X-ray and optically selected AGNs, Ross et al. for SDSS QSOs, and Groom et al for 2dF QSOs are plotted with markers specified in the legend in this 
figure. The horizontal error bars represent the redshift range of the samples. 
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14 
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_ 10 
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6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 

log(MBH/Msoi) 

Fig. 12. Cross-correlation lengths derived for each redshift matched samples are plotted as a function BH mass. The data point of the same redshift groups 
are connected with a line. Filled circles are for zl redshift group, filled triangles and boxes are for z 2 and 23 redshift groups, respectively. 


redshift matched sample 



15 


10 


5 


0 

0.0 0.2 0.4 0.6 0.8 1.0 

Redshift 

Fig. 13. Cross-correlation lengths derived for each BH mass matched samples are plotted as a function of redshift. The data point of the same mass groups 
are connected with a line. Filled boxes are for M75 mass group, filled circles and triangles are for M82 and M90 mass groups, respectively. 




























20 


Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0 




I^opt-iR l^opt-IR l^ppt-IR 

Fig. 14. Distributions of Dopt-m parameter for each redshift and mass groups. Fiiied circies represent observed data points, soiid lines are fitting result of 
the two component model, red solid and blue dashed lines are for red and blue galaxy component, respectively. 



<Log(MB„/M„,)> 



Fig. 15. Left panel: Red galaxy fractions derived from the two component fitting shown in Figure 14 are plotted as a function of BFI mass. The results of the 
same redshift groups are connected with a line. Circles are results for redshift z1 group, and triangles and boxes are for redshift z2 and z3 groups, respectively. 
Right panel: Normalized excess density derived from the two component fitting. Filled and open circles are results for blue and red galaxies respectively for z1 
redshift group. Filled and open triangles are results for z2 redshift groups with the same blue/red galaxy assignment. Filled and open boxes are results for z3 
redshift groups with the same assignment, and the arrows represents the upper limit for blue galaxies at redshift z3. 
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Fig. 16. Distributions of Dopt-m parameter for redshift matched samples. Filled circles represent observed data points, solid lines are fitting result of the two 
component model, red solid and blue dashed lines are for red and blue galaxy component, respectively. 
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Fig. 17. Left panel: Red galaxy fractions derived from the two component fitting shown in Figure 16 are piotted as a function of BH mass. The results of the 
same redshift groups are connected with a line. Circles are results for redshift z1 group, and triangles and boxes are for redshift z2 and z3 groups, respectively. 
Right panel: Normalized excess density derived from the two component fitting. Filled and open circles are results for blue and red galaxies respectively for z1 
redshift group. Filled and open triangles are results for z2 redshift groups with the same blue/red galaxy assignment. Filled and open boxes are results for z3 
redshift groups with the same assignment. 






Fig. 18. Comparison of absolute magnitude distributions between higher (red histogram) and lower (blue histogram) mass groups. The ratios of high/low are 
plotted in the panel below the histograms. Dashed lines shown in the ratio plot is a result of a linear fit in linear-log scale. 





















































Excess density ratio (high mass / low mass) 
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Fig. 19. Comparisons of the distributions of average number densities of gaiaxies at AGN redshift between higher (red histogram) and lower (blue histogram) 
mass groups. The left (right) panel is the comparison at redshift z1 (z2+z3). 



I^IR 


Fig. 20. Comparison of the ratio of absolute magnitude distributions. Filled blue circles represents the ratio of M82 and M65+M75 at redshift z1. Open red 
circles represents the ratio of M90 and M75+M82 at redshift z2 and z3. Fitting results with a power law function are also shown with the solid lines. 




































































